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ABSTRACT 

We present first results from a large-area (~80°x20°), sparsely sampled two- filter (B,R) imaging 
survey towards the Canis Major stellar over-density, claimed to be a disrupting Milky Way satellite 
galaxy. Utilizing stellar colour-magnitude diagrams reaching to B^22 mag, we provide a first delin- 
eation of its surface density distribution using main sequence stars. It is located below the Galactic 
mid-plane, and can be discerned to at least b=-15°. Its projected shape is highly elongated, nearly 
parallel to the Galactic plane, with an axis ratio of > 5:1, substantially more so than what Martin 
et al. originally found. We also provide a first map of a prominent over-density of blue, presumably 
younger main sequence stars, which appears to have a maximum near [l,b ~ 240°, -7° ] and extends in 
latitude to b~ -10°. The young population is markedly more localized. We estimate an upper limit on 
the line-of-sight depth, tri os , of the old population based on the main sequence width, by comparing 
with a simple stellar population (Pal 12), obtaining <ri os < 1.8 ± 0.3 kpc, at an adopted D = 7.5 ± 
1 kpc. For the young stellar population, we find <j\ os < 1.5 kpc. The overall picture presented is one of 
a young stellar population that is less extended, both in terms its line-of-sight depth and angular size, 
than the older population. In the literature three explanations for the CMa over-density have been 
invoked, namely (a) a partially disrupting dwarf galaxy on a low-latitude orbit, (b) a projection of a 
warped outer Galactic disk, and (c) a projection of an out-of-plane spiral arm. While the data provide 
no firm arguments against the less well-defined third scenario, they have clear implications for each of 
the others: (a) We infer from the strong elongation of the over-density in longitude, and simulations 
in the literature, that the CMa over-density is unlikely to be a gravitationally bound system at the 
present epoch, but may well be just a recently disrupted satellite remnant. A complexity for the 
satellite origin may arise from the 'flattening' of the young MS population, which is possibly more 
pronounced than the older one. (b) Based on modeling, we find that the line-of-sight depth of the MS 
over-density in old stars is clearly inconsistent with published locally axi-symmetric descriptions of 
the warped Galactic disk, such as those considered in Momany et al. (2006). Without detailed mod- 
eling, the data-set itself does not provide sufficient leverage to distinguish between an interpretation 
as sub-structure in the warped outer Galactic disk or a disrupted satellite. 

Subject headings: galaxies: dwarf — galaxies: interactions - galaxies: individual (Canis Major) - 
Galaxy: evolution - Galaxy: stellar content — Galaxy: structure 



1. INTRODUCTION 

The current paradigm in which large disk galaxies like 
our Milky Way form (e.g. White & Rees 1978; White 
& Frenk 1991), is based on the successive coalescence 
and accretion of smaller systems of dark matter, stars, 
and diffuse interstellar gas into larger assemblies. Some 
of these smaller systems may be satellite dwarf galaxies 
that can be disrupted by tidal shocks and evaporation 
in a Galactic potential, and can spawn tidal tails. Nu- 
merical simulations exist that model distinct tidal stel- 
lar streams in and around large galaxies, indicating that 
at r> lOkpc such streams should remain detectable as 
coherent stellar over-densities for billions of years (John- 
ston et al. 1999; Ibata & Lewis 1998; Martinez-Delgado 
et al. 2004; Penarrubia et al. 2006). 

Recently, much work has focused on a possible low- 
latitude stellar stream, the so-called Monoceros stream, 
encircling the Milky Way (MW) at a galacto-centric dis- 
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tance of around 20 kpc, which was found through colour- 
magnitude diagrams (CMDs) from SDSS 1 (Newberg et 
al. 2002; Yanny et al. 2003); also see Ibata et al. (2003), 
Conn et al. (2005) & Martin et al. (2006). For stellar 
over-densities at low galactic latitudes it is not obvious 
a priori whether such a stream is of external origin, i.e. 
the tidal debris of a now (partially) disrupted satellite or, 
alternatively, a distorted part of the pre-existing outer 
stellar disk. 

Arguments in favour of the dwarf satellite hypothe- 
sis have come from dynamical modeling (Penarrubia et 
al. 2005), showing that it is possible to find a plausible 
model of the Monoceros stream that explains all detected 
parts of this low-latitude stream as the wrapped tidal de- 
bris tails of a disrupting (model) dwarf galaxy. In this 
dynamical model the position of the CMa over-density 
is not an initial constraint, and yet the main body of 
the disrupting model satellite happens to located in the 

1 Sloan Digital Sky Survey 
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direction of CMa at the present epoch. All kinemat- 
ics, including the subsequently determined proper mo- 
tion (Dinescu ct al. 2005) of main sequence (MS) stars 
in a small area (0.25 deg 2 ) towards the CMa stellar over- 
density are consistent with the dynamical model. Addi- 
tional empirical support appears to come from the ap- 
parent bifurcation of suspected Monoceros stream stars 
at right ascension ~ 120° in Fig.l of Belokurov et al. 
(2006), as could be expected in general based on the 
above-mentioned Penarrubia et al. model. Other ev- 
idence comes from the discovery of an over-density of 
possible Monoceros stream RR Lyrae stars (Vivas et al. 
2006), whose distances closely match the model expecta- 
tion in Penarrubia et al. (2005). 

The position and survival of the progenitor of this stel- 
lar stream is still under debate. The best candidate 
is a seemingly well-defined stellar over-density of stars 
discovered in the direction of Canis Major (CMa) us- 
ing 2MASS red giants (Martin et al. 2004a). This dis- 
covery has spawned a lively debate in the literature on 
whether this apparent over-density of red giants could 
be part of the distorted outer stellar disk (Momany et 
al. 2004; 2006; see also Rocha-Pinto et al. 2006) or an 
accreted and possibly disrupting satellite [Martin et al. 
2004b; Martinez-Delgado et al 2005 (MD05); Bellazzini 
et al. 2006 (B06)]. However, a firm large-area kinemat- 
ical, spatial and chemical link between the three known 
stellar components of the overall debate (i.e. the young 
and old stellar over-densities, and the Monoceros stream) 
is lacking. 

From deep, visible-band photometry there are signs of 
at least two different star formation episodes in the direc- 
tion of CMa [Bellazzini et al. (2005), B05; MD05; Car- 
raro et al. (2005)]. For the sake of clarity, we illustrate 
this in an annotated colour-magnitude diagram, Fig. ^ 
which corresponds to a field near the presumed centre 
of the CMa over-density, and where the data are taken 
from this survey (see Sec. EJ- Whether the most recent 
burst of star formation activity occurred 1-2 Gyr ago, as 
touted in B05, has been debated, in favour of a much 
younger stellar population of < lOOMyr-old stars (Car- 
raro et al. 2005). The CMa over-density itself, which is 
taken to lie near D Q =7.5kpc (B05,MD05), comprises a 
predominantly older population of stars (4-10 Gyr; B05) 
and Moitinho et al. (2006) suggest that it may be con- 
sequence of viewing a local arm-like MW stellar sub- 
structure in projection. For a different and independent 
study of the stellar populations, we will make use of de- 
tailed CMD-fitting in a forth-coming paper (de Jong et 
al., in prep.). 

To determine the full angular extent of the young 
and old main sequence stars towards CMa and to test 
whether they stem from a satellite galaxy, we conducted 
a large-area (^80°x20°), sparsely sampled visible-band 
imaging survey of the CMa region, drawing on the wide 
field imager (WFI) at the MPG/ESO 2.2 m telescope on 
La Silla. The survey sub-area considered in this paper is 
230° < 1 < 260°, -20° < b < 15°. 

1.1. Aim and design of the WFI survey 

The principal goal of our imaging survey is to provide a 
database that can address ultimately whether the CMa 
stellar over-density results from a dynamical distortion 
of the outer Galactic disk or whether it is the remains 



of a formerly large dwarf galaxy. As a step towards this 
goal, we aim to map its full angular extent through its 
star-count profile. 

Recent reports (B06) based on 2MASS red clump stars 
suggest that the discernable part of the CMa over-density 
may extend across several lOOdeg 2 on the sky. Our 
survey goal is to cover this over-density through sparse 
sampling, but with relatively deep visible-band imaging. 
Specifically, we want to obtain CMDs that reach well 
below the turn-off magnitude of an ancient stellar popu- 
lation; our target limiting magnitude is B, R = 23 mag 
(S/N = 10). We chose a wide colour baseline (B-R) to 
have sensitivity to age- and/ or [Fe/H]-dependent CMD 
features (young MS, old MS turn-off and the red clump). 
The key difference (and merit) of the present CMa survey 
over previous large-area surveys, which have been based 
exclusively on 2MASS red-giant stars and MW models, is 
that we base our analysis on colour-magnitude diagrams 
that exhibit a well-defined stellar main sequence at a de- 
fined distance range. This provides an independent and 
high-contrast way of tracing the full angular extent of 
the CMa over-density. The chief disadvantage of visible- 
band data is that Galactic extinction in many parts of 
this area is high, Ay >0. 5 mag. 

1.2. Aim of this paper 

In this paper we present the data from the first ob- 
serving phase of our imaging survey, covering the region 
around [l,b = 240°, -8°] that had originally been iden- 
tified as the over-density in Martin et al. (2004a). As 
it turns-out, these data provide only partial spatial cov- 
erage of the over-density. Therefore, we refrain from a 
wide-spread quantitative comparison with Galaxy star- 
count models, except (i) to use a Besancon model CMD 
as a visual guide for the reader to help interpret the gen- 
eral morphology in our control CMDs, and (ii) to aid the 
discussion in Sec. 17.21 

We report our observations, and present the data and 
their reduction in Sec. El The issue of photometry com- 
pleteness and uncertainties is assessed in Sec. 12.11 The 
crucial issue of reddening by dust extinction and its cor- 
rection is presented in Sec. 12.21 In Sec. El we explain the 
choice of control fields, which is followed by a descrip- 
tion of the stellar content and morphology of the CMDs 
in Sec. 01 In Sec. El we provide a star-count analysis of 
young and old CMa main sequence stars. We estimate 
an upper limit on the line-of-sight size of the CMa stellar 
overdensity in Sec. El We briefly discuss our results in 
Sec. Hand summarize the key results in Sec. El 

2. OBSERVATIONS, DATA AND DATA REDUCTION 

All observations were carried out in service mode with 
the wide-field imager (WFI) on the 2.2-m ESO/MPG 
telescope at the La Silla observatory (Chile) from De- 
cember 9 to 20, 2004. The WFI's field-of-view covers 
0.25deg 2 , sampled at 0.238" pixel -1 . Single exposure 
images were taken in the B- and R-bands, at 100 s per 
pointing and filter. A total of 59 pointings (Table Q 
provides a sparse map of the stellar over-density region 
originally identified in Martin et al. (2004a). B- and 
R-band pointings typically differ by less than 2.5arcsec, 
which leads to negligible field-to-field differences in the 
(effective) imaging area. Overscan, bias, flat-field cor- 
rections, and an astrometric correction were performed 
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using a pre-reduction pipeline (Schirmer et al. 2003). 
We obtained stellar photometry using DAOPHOT tools, 
available in IRAF. 2 For each WFI pointing, we employed 
point spread function (PSF)-fitting using a spatially vari- 
able (order 2) moffat function of exponent 2.5, derived 
from 100 stars. The DAOPHOT/ ALLSTAR PSF-fitting 
task is run on each frame separately. As we have sin- 
gle exposures per pointing, cosmic rays are effectively 
filtered out during the PSF-fitting and the matching of 
the B- and R-band photometry lists. Objects in each list 
are then assigned Galactic coordinates and are matched 
using a separation tolerance of 1.7arcsec. This large tol- 
erance ensures that the faintest stars, which have the 
poorest positional errors, are matched; and is tolerable 
owing to the sparseness of the fields surveyed. 

For the photometric zero-points, we firstly tie the pho- 
tometry from each observing night to a common inter- 
nal (i.e. instrumental) zero point using the zero-point 
of each PSF (computed from its flux and magnitude). 
This is followed by the transformation of the instrumen- 
tal magnitudes into the standard Johnson-Cousins pho- 
tometric system based on observations of either standard 
stars from Landolt (1992) [for Dec. 17, 18] or secondary 
standards from Galadi-Enriquez et al. (2000) [on Dec. 
20.] or both [Dec. 9], taken before or after our obser- 
vations. For all standard star observations, atmospheric 
conditions were photometric. The calibration photome- 
try, obtained from large (r=2.1 arcsec) circular apertures, 
is used to transform our instrumental magnitudes (b, r) 
into the Johnson-Cousins standard system (B,R). We es- 
timate the colour terms (see Eqiis-^and 0) using 20 cal- 
ibration stars from Dec. 9 with standard colours (B-R) in 
the range ~ 0.8 to 1.6 mag. While we do not have enough 
calibration stars to determine accurate colour terms for 
each night separately, CMa observations were always per- 
formed under photometrically clear conditions, and so we 
use the colour terms from Dec. 9 for each of the observa- 
tions. For a determination of the photometric zero-point 
offset (s), we have 4 to 20 calibration stars each night, 
except for Dec. 15 & 16., for which the calibration from 
Dec. 18 is used. With the instrumental fluxes already 
scaled to photons per second, the transformation equa- 
tions for Dec. 18, for example, are 

B = -b + 24.86(±0.01) + 0.17(±0.01) (B - R) (1) 

R= -r + 24.48(±0.04) - 0.07(±0.03) (B - R) (2) 

where the instrumental magnitudes (b,r) are corrected 
for atmospheric extinction. Wc adopted the atmospheric 
extinction co-efficients given in the ESO/WFI web-pages; 
those are 0.24 and 0.09 for the B- and R-bands respec- 
tively. The error values given in Eqns.nand [21 are the 
rms uncertainties. 

The survey photometry is filtered using quality-control 
parameters computed by ALLSTAR in DAOPHOT, re- 
taining only stars with acceptable CHI and SHARP 
parameters as well as formal magnitude errors ctb < 
0.15 mag and ctr < 0.15 mag. Galactic coordinates and 
the number of selected stars per field is given in Table 

2.1. Completeness and photometry errors 

To assess the point-source completeness, we performed 
artificial star tests for a representative set of six fields. 

2 Image Reduction and Analysis Facility. 



Specifically, we added a total of several thousand stars, 
in groups of 1300, in the (instrumental) magnitude range 
15.5 to 25.5 mag, to each image. The fake stars that we 
added at random locations to these largely uncrowded 
fields increase the number-density by up to 10%, hence 
barely altering crowding effects. Stars injected in the 
B- and R-band images did not have the same locations. 
We then determined the completeness fractions, defined 
as the ratio of recovered artificial stars to the number 
of the injected ones, at 0.75 mag-wide intervals. We 
brighten the B- and R- fake star magnitudes by 0.42 and 
0.13 mag, respectively, to shift them approximately from 
instrumental magnitudes to the flux-calibrated magni- 
tude scale. Also, as we use de-reddened CMDs in this 
paper, we shift the the magnitude scale of the complete- 
ness curves accordingly and denote the shifted B- and 
R-band magnitudes by B' and R' respectively in Fig. [21 
The completeness at (l,b) = (241.5°,-6°) is above 90% 
for B' or R < 20.6 mag (see Fig. EJ. In the matched 
B- and R-band photometry lists the completeness is less, 
and is typically above ~ 80% at B' < 20.6 mag. This 
80% completeness limit varies by ± ~ 0.5 mag across the 
area surveyed at |b|>6°. 

The completeness at each pointing is reduced by gaps 
between detector chips, bad columns and pixels that are 
the same for each field. Completeness at the faint magni- 
tude end varies from field-to-field for a number of reasons: 
differences in the seeing, stellar crowding, the cumulative 
effective of saturated stars (together with their stray light 
and ghost images, caused by internal reflections), cosmic 
rays, bright galaxies, and bright, nearby moving objects. 
Based on the database of artificial star photometry, we 
can assess the photometry errors (e.g. see Fig.[3J). For ex- 
ample, 68% of all artificial stars and, by inference, all de- 
tected stars with B'~22.6 and R'~22mag, at |b|> 6deg, 
have a photometry error typically below 0.05 mag, and 
smaller still at brighter magnitudes. 

2.2. Extinction and extinction correction 

Since all our target fields are at low Galactic latitudes, 
the line-of-sight dust extinction (and, to a lesser extent, 
its dispersion) within each field is of great importance to 
the analysis in this paper. Schlegel, Finkbeiner & Davics 
(1998; SFD98) dust maps provide extinction estimates 
for high latitude fields, but are likely to be inaccurate 
near the Galactic mid-plane, particularly at |b|< 10°, 
according to SFD98. 

As we cannot rely on SFD98 corrections alone, we per- 
form an empirical foreground correction of the CMDs in 
two steps. For the first step we take the SFD98 dust 
maps, interpolated on a star-by-star basis. Adopting a 
standard Galactic extinction law (Cardelli, Clayton & 
Mathis 1988) we then determine modified de-reddening 
values based on the prescription given by Bonifacio, 
Monai & Beers (2000; their Eqn. 1), and apply those 
to the photometry; the resultant differential extinction 
data is denoted simply as E(B-V)sfd* hereafter for the 
sake of clarity. However, the colour of the resulting CMD 
morphology in different fields shows that this cannot be 
the full reddening correction. Therefore, we apply a sec- 
ond correction step, enabled by the prominent main se- 
quence feature in the resulting CMDs and based on the 
assumption that the metallicity and age of the stellar 
population that causes this CMa MS over-density (e.g. 
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see MD05) does not vary much from field-to-field (though 
its distance may vary). With this assumption, the de- 
reddened colour of the main sequence turn-off (MSTO) 
stars should be the same in each field. Therefore, we 
shift the stars in each CMD along the reddening vec- 
tor to match the colour of the MSTO region 3 at [l,b = 
242.5°, -9.0°; E(B-V) S fd» =0.148 mag], which is a rela- 
tively low extinction field with a significant CMa signal. 
This shift is of course a shift both in colour and magi- 
tude. To obtain each MSTO colour we make a coarse es- 
timate (number- weighted centroid) at Bo = 18.5-19 mag, 
and then re-compute it in a 0.5 mag- wide colour interval, 
centred on the coarse estimate, binned at 0.03 mag in- 
tervals. Outlying bins are ignored by setting a bin-count 
threshold, taken as the median bin-count. When testing 
the de-reddening procedure, we identified obvious fail- 
ures through a visual inspection of the CMDs, focussing 
on the colour of the MSTO region. As the CMa MSTO 
is not well defined in every field, we do not apply this ad- 
ditional de-reddening step if E(B-V)sfd* < 0.11 mag or 
if the field is outside -14.7° < b < 4°. De-reddened fields 
near b~-15° and also at b > 4° have no significant CMa 
content, and therefore can appear too red by < 0.1 mag. 
However, as the density of stars near (B-R)o =1 mag ap- 
pears to vary smoothly in our CMDs at such Galactic 
latitudes, this effect is only a minor concern throughout 
this paper. Magnitudes and colours de-reddened in the 
above-mentioned fashion are denoted by B and (B-R) 
respectively in this paper, and by Bo and (B-R)o other- 
wise. 

The (MSTO) colour shifts resulting from the de- 
reddening process, detailed above, are also given in Ta- 
ble n those values tend to grow with increasing Galactic 
longitude, and tend to diminish away from the Galac- 
tic mid-plane. The values suggest that the SFD98 data 
both over- and under-estimates the foreground extinction 
across our survey area. In order to delimit the impact of 
de-reddening errors further, we take an (ad hoc) extinc- 
tion threshold E(B-V) SFD *=0.30mag. E(B-V) SF d* data 
is given in Tabled for each field. 

A final issue is the reddening dispersion in each field. 
The only attempt at compensating for it in this paper 
is through the first de-reddening step, detailed above, 
which we apply on a star-by-star basis. As the actual 
reddening dispersion is unknown, we can only proceed 
on the assumption that the predominant variation from 
field-to-field is in the (median) foreground dust extinc- 
tion, which we have attempted to correct for, and that 
there are no significant differences in the residual (fore- 
ground) dispersion from field-to-field. 

3. CONTROL FIELDS 

The overall density and the line-of-sight distribution of 
the outer Galaxy's stellar constituents (halo, thick/thin 
disk) vary with position on the sky. As we do not pre- 
cisely model how the known (locally) axisymmetric con- 
tributions vary in the area of sky sampled in this pa- 
per, we try to minimize the dependence of our results 
on the current generation of synthetic MW models. We 
therefore adopt an empirical estimate of the MW com- 
ponents' contribution to our CMDs. To illustrate how 

3 We refer to the blue near- vertical/arc-like over-dense region at 
B ~18.5-to-20 mag (E.g., see Fig|T) 



we select such 'control fields', we plot in Fig. 0a sub-set 
of CMDs for three fields of similar |b|, from above and 
below the Galactic mid-plane. It is immediately appar- 
ent from this that the CMa over-density is not present 
in each field, and those without obvious visible evidence 
of the CMa-over-density are termed control- fields. As 
we have control fields both above and below the Galac- 
tic mid-plane, we take the simple approach of using the 
field at (l,b) =(240°, +8°) for the latitude range -14° < 
b < 8°, and take the (240°,-20°) field when b <-14°. At 
b> 8°, the field at (240°, +15°) is the control field. We 
stress that the full imaging survey will provide a larger 
array of control fields for a better assessment and usage 
of them in future analyses. 

4. COLOUR-MAGNITUDE DIAGRAM MORPHOLOGY AND 
THE STELLAR POPULATIONS TOWARDS CMA 

4.1. CMD morphology and stellar content of a control 

field 

Fig. 01 shows the CMD of a field that we take as a con- 
trol field from the current phase of the imaging survey, 
together with the Besancon MW model counterpart 4 . It 
was chosen as a control-like field because it shows no 
obvious presence of the old CMa MS reported in BOS 
and MD05, based on a visual inspection. The Besan- 
con model is meant to be a description of a MW galaxy 
without inhomogeneities (e.g. spiral arms). As such it 
is very useful for exploring what types of stars at what 
distances and from which MW components (halo, disk, 
spiral arms) contribute to different parts of the CMD. 

The CMD in Fig. ^comprises a range of stellar popula- 
tions at widely differing distances. The most prominent 
feature of the control field and its model counterpart is 
the so-called blue-edge of MSTO stars. The colour of this 
edge at B < 19 mag is influenced by metal-rich thick disk 
stars and probable metal-poor halo MSTO stars. It is ~ 
0.1 mag bluer at fainter magnitudes, where it is domi- 
nated by metal-poor, outer-halo MSTO stars (spectral 
type F). We note (1) that apart from possible photome- 
try scatter, some of the stars at (B-R)o <-0.4 and Bo> 
19 mag may be nearby (<1.5 kpc) white dwarf stars. Blue 
horizontal branch stars in the halo and thick disk would 
also have such blue colours. (2) At Bo < 15 and (B-R) 
< 0.6 mag in the Galactic model CMD, we see that one 
can expect nearby thin/thick disk MS stars within about 
2 kpc. Further, (3) there is the generally truncated ap- 
pearance of the stars distributed on the red side of the 
CMD, including a red plume of late-type (M) MS stars, 
which can appear tilted in some fields. As can be seen 
in Fig. EJbottom) , those stars are nearby thin and thick 
MS stars within about 2 kpc, and to a lesser extent, red 
giant branch stars farther out. 

Fig. El illustrates that the structure of the CMD in 
the Galactic North and in the South are different over 
the full latitude range at 1~ 240°. A comparison with 
Fig. Ufbottom) shows that the structure of the control 
field CMDs can be reproduced using a standard Galactic 
model such as the Besancon model, but below the mid- 
plane an extra population is needed to create the CMa 
over-density, which is attributed to the MS of a distinct 
stellar system near D =7.5kpc in B05 and MD05 (their 

4 Default Besancon model via the online web-page at bison. obs- 
besancon.fr/modele/; Robin et al. (2003) 
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Fig.l); see Sec. O also. 

4.2. CMD morphology and stellar content of CMa fields 

Fig. ^ explains what one expects to see in a typ- 
ical CMa CMD. Comparing the control fields (i.e. 
b=+15,+8,-20° in Fig.© and this CMa field, we can now 
discern more robustly the same CMa features picked out 
by MD05 for a pointing near the presumed centre of the 
CMa over-density. A prominent feature is the old MS 
(B05; MD05), which has a blue, arc- like (possible MS 
turn-off) region near B = 19 mag. At brighter magni- 
tudes there is a blue plume of young MS stars at [B ,(B- 
R) ] ~ [-0.2, 15] to [0.5, 18] mag. Further, there is the 
hint of a red-giant branch in many CMDs, that is visible 
from [B ,(B-R) ] - [15.5, 1.5] to [17.5, 1.2], and it may 
contain some CMa stars. 

Fig. El to 0show a large representative set of 15 tiles 
from our sparse map and provide a view of the stellar 
content in the CMa stellar over-density. These figures 
cover a good range in galactic longitude and latitude, 
and provide a first qualitative picture of the stellar con- 
tent of the CMa stellar over-density over a large angular 
area. We see that the old MS is a very high-contrast 
feature in several of those tiles, e.g. Fig. ^centre panel) 
or Fig. [7| (top-right panel) . The density of old MS stars 
exhibits a prominent variation with latitude, decreasing 
away from the Galactic mid-plane (Fig. EJ). Comparing 
with CMDs along the longitude direction (Fig[7J), a dis- 
tinct elongation in the old MS stellar population is imme- 
diately apparent. In contrast, we see that the young MS 
population is less evident and less populous at greater 
longitudes. From Fig. El and Fig. [7| it is apparent that 
the old MS is still present at b=-15°, while the younger 
population is still present at b=-10°. 

The number-density and shape of the young MS star 
population varies among the CMDs. Fig. El shows that 
there are dense plumes with a broad B-band spread, e.g. 
[l,b = 240.6°, -6.9°], but there are also narrow, dense 
plumes, e.g. [l,b = 237.4°, -8.1°] (see Fig. 01. The di- 
versity in the shape (B-band width and density) of the 
young MS is reminiscent of the Small Magellanic Cloud 
(e.g. see CMDs in Noel et al. 2005), where there has 
been on-going and bursty star formation over the past 
few billon years (Harris & Zaritsky 2004). The blue MS 
in the CMa CMDs may contain some blue straggler stars. 
It is likely to be dominated by young MS stars, which is 
how we refer to this population for the remainder of this 
paper. Carraro et al. (2005) report that this young MS 
is also observed in all of their fields, closer to the Galactic 
mid-plane, and that this suggests it is associated with the 
MW spiral galaxy. However, this interpretation may be 
complicated by the observed differences in its morphol- 
ogy (density and magnitude width) at different galactic 
longitudes (Fig. CJ). 

This young (CMa) MS star plume should not be con- 
fused with the other similarly blue swath of foreground 
(< 2 kpc) white dwarfs at fainter magnitudes in most of 
our CMDs. Exceptions occur close to the Galactic mid- 
plane, owing to the possibly inaccurate de-reddening of 
such nearby stars. The young (CMa) MS stars are also 
different from young (<2-3Gyr) foreground (< 4 kpc) 
thick/thin disk MS stars that can occur at similar and 
brighter magnitudes than the young MS stars [e.g., at 
Bq < 19, (B-R) < 0.4 mag] marked in Fig-ffl Such MW 



disk stars can be seen at b=5.0° and b =-5.9° (1=240°) 
in Fig. EJ for example. Based on the model CMD [see 
Fig. E£ top)], those stars are foreground thin/thick disk 
MS stars at D <2kpc. That such MS stars are not de- 
tected in every field may be a consequence of field-to-field 
differences in the saturation magnitude, which in turn 
depends on foreground dust extinction, 5 and observing 
(atmospheric) conditions. 

5. MAPPING CMA ON THE SKY 

5.1. Method - old and young MS 

To estimate the old MS star-counts that are at- 
tributable to the CMa over-density at each pointing, we 
apply a simple analysis of the de-reddened photometry 
using CMD extraction boxes. We estimate the number 
of old CMa MS stars, NcMaMS) m the relevant extraction 
box (see Fig. EJ at a given pointing (l,b) as follows: 

AT r A T -^control field;MSbox Ar -i 

JVCMaMS = 7 [ iv MS box — "TT ^V re f.boxJ, 

-^control ficld;rcf .box 

(3) 

where N con t ro i fi e id;MSbox is the star-count in the MS 
extraction box in the relevant control-field. Nms box and 
Nrof.box are the star-counts at a given (l,b) in the CMa 
old MS box and a reference box respectively. These boxes 
are outlined in Fig-IHl 7 is the scale factor needed to com- 
pensate for the fraction of area (actually star-counts) at 
16 < B < 20 above the adopted extinction threshold, 
E(B-V)sfd* =0.30 mag, in each field. 7 = 1 for the ma- 
jority (> 90%) of these fields, and l-to-2 otherwise. A 
density estimate and its (random) uncertainty for each 
pointing is then obtained through data re-sampling. I.e., 
100% of the stars are selected at random, allowing re- 
peated selection, and a MS width estimate is recorded. 
After repeating this 100 times, we fit a Gaussian function 
to a histogram of the density estimates. The mean value 
from the fit is taken as the density estimate and the width 
parameter from the fit is taken as an estimate of the as- 
sociated (random) error. Both of them are recorded in 
Table [U 

We also map the surface density of young MS stars. 
We do this by counting stars in the extraction box over- 
laid on the young MS in Fig. |H1 While there are many 
fewer young MS stars than old MS stars, any contamina- 
tion of the former by field stars appears to be negligible, 
based on either a visual inspection of our CMDs or the 
Besancon model CMD in the same direction. For that 
reason we do not subtract a background estimate. Each 
density estimate and its error (Table is determined 
through data re-sampling. 

5.2. Results: Comparison of young and old MS 
star- count profiles 

Fig. El shows the estimated number of young and old 
MS stars per WFI field against Galactic latitude and 
longitude. These are the key results: 

The old MS stellar over-density is highly elongated in 
Galactic longitude in the area of sky surveyed [also see 
Fig. El (bottom)]. Although complicated by low-latitude 
extinction (and missing sky coverage) , the position of its 

Sufficient foreground extinction can cause otherwise unde- 
tected (saturated) stars to be detected by dimming them. 
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peak density appears to be located at (1, b) = (>240°, 
>-7°), based on a visual inspection. In turn, this sug- 
gests that CMa lies predominantly (e.g. >60%) below 
the Galactic mid-plane, as found in the 2MASS red-giant 
star analysis in Martin et al. (2004a). The survey data 
constrain the projected aspect ratio of the old stellar 
over-density. There is little surface density gradient in 
longitude across the survey area and it is therefore quite 
reasonable to infer that its FWHMi exceeds the survey 
width, i.e. > 27°. In the latitude direction, the angular 
width estimate of the CMa density profile is complicated 
by the fact that the density maximum can only be brack- 
eted to be at b= < -7° in Fig. El Taking the turn-over 
to be at b~-7°, the FWHM b is - 6°, based on a visual 
inspection, and CMa's projected aspect ratio is therefore 
> 5:1. For the young MS, taking FWHM, > 20°, con- 
servatively, and FWHMb ~ 2°, would give a projected 
aspect ratio of > 10:1. 

The young and old MS star populations overlap on 
the sky, but the young stars are markedly more localized 
in latitude. The detection of a compact distribution of 
young stars in the longitude direction is especially signif- 
icant because of the essentially uncontaminated sample. 

The over-density of young MS stars is not as extended 
in projection as that of the older stellar population [see 
Fig- EH- but is possibly more 'flattened.' The young MS 
stars exhibit a maximum density at (1, b) ~ (240°, -7°), 
with a drop-off in their counts at 1<240° and 1> 240°, as 
can be verified qualitatively from a visual inspection of 
the CMDs in Fig. [3 The same basic conclusions regard- 
ing those stars can also be drawn from Fig. which 
shows a sky-map of their number-density. Lastly, there 
are wiggles in the young and old MS profiles, but there 
is the possibility that they result from having a small set 
of control fields, and/or unknown reddening dispersion. 

6. LINE-OF-SIGHT DEPTH OF CMA AT (L,B) = (242.5°, -9°) 

MD05 argued that the quite narrow distance range 
of the CMa stars points towards a disrupting satellite, 
rather than a flaring or warping of the outer Galactic 
disk. Here, we explore the line-of-sight extent of CMa 
by estimating the MS width in the magnitude direction 
for (l,b) = (242.5°,-9°), a low-extinction field near the 
presumed centre of CMa. As this estimate is expected to 
be the same over the CMa over-density we compare with 
the estimated MS width in MD05 for (l,b) = (240° ,-8°). 

6.1. Estimating the MS width 

We model the B-band star-count distribution, f tota i 
(Bq), in a given CMD colour slice as the linear sum of 
two components, namely a smooth (underlying) distri- 
bution of stars from the Galactic stellar halo, thick and 
thin disks, fbkg(B ), and the CMa MS, fMs(B ): 

ftotai (Bq) = afbkg(Bg) + AfMs(Bo), (4) 

_where i M s(B ) - [l./(x/2^7 Bi )] EXP[-(B - 

B ) 2 /(2erg, )], and the coefficients (a, A) are scale 

factors. We describe the background distribution of 
stars, fbkg(Bg), as a second order polynomial, as detailed 
below, and allow only its amplitude to vary. 

In order to have a compact B-band distribution of CMa 
stars, while having a wide colour slice, 1.0< (B-R)q < 



1.1 , for signal-to- noise reasons, we consider the distri- 
bution (B^ „ = ) B' Q -3.5[(B-R)(, - 1.05], instead of B' 
alone. The slope is chosen to closely match the slope of 
the old MS in our CMDs. The method is illustrated in 
Fig. HU for (l,b) = (242.5°, -9.°), a field near the pre- 
sumed centre of CMa, and the reference field (l,b) = 
(240°, +8°). The analysis differs from MD05, who had 
no control field, and a fainter magnitude limit, which al- 
lowed them to consider a narrower colour slice further 
down along the MS. The functional form of fi,k g is given 
in Fig. ^2 (top- left panel), with its amplitude scaled to 
match fjvis at Bq „ = 18-19 mag, binned at 0.25 mag in- 
tervals. After subtracting fbkg, h/is is fitted. The er- 
ror introduced by the binning and shot noise is assessed 
through data re-sampling. For the observed MS width, 
we obtain <<tb' > = 0.49 ±0.05 mag. 

In the above construction, we have assumed a sym- 
metrical distribution of star-counts N(Bq„), which ap- 
pears to be a reasonable first approximation based on 
a visual inspection of Fig. II If bottom-left ). To convert 
the estimated width of the observed MS, < a' Bo >, into 
a physical line-of-sight (l.o.s.) depth, we firstly sub- 
tract in quadrature the contribution from the photo- 
metric uncertainty, and the intrinsic width, based on an 
(assumed) stellar population with a zero distance range, 
thereby obtaining a m . In principle, a reddening disper- 
sion should also be subtracted, but our only available 
way of attempting this is already implemented through 
the star-by-star de-reddening detailed in Sec. 12.21 We 
then transfer a m to a physical depth estimate by making 
use of m - M = -5 + 5Logio D©, where D© (in par- 
sees) is the heliocentric distance of a star of magnitude 
m. The determination of a\ os is complicated by the fact 
that the adopted fMS m ay not be completely accurate, 
as remarked in Lopez-Corredoira (2006). A good first 
approximation for the physical depth, in terms of the 
line-of-sight la 'radius' or depth, is this: 

a los « Dq [10 ct »/ 5 - 1] (5) 

where D Q is the distance to the barycentre of the CMa 
over-density along the line-of-sight. It is taken to be 
7.5±lkpc, a mid-way estimate between values in the lit- 
erature [7.2±lkpc, B06; 5-8 kpc, MD05]. 

6.2. An upper limit on the line-of-sight (l.o.s.) depth of 

CMa 

To estimate the depth of the old and young CMa over- 
densities, we compare the width of its MS to that of a 
reference population. We now consider the old MS pop- 
ulation and return to the young MS later below. For the 
old over-density, the reference is a simple stellar popu- 
lation with essentially no depth (i.e. cri OSirc f/D0 iro f << 
1). Thus, by construction, we do not account for a likely 
spread in the chemical composition of the (old) CMa stel- 
lar population and can provide therefore only an upper 
limit estimate of its line-of-sight extent. 

For the old reference stellar population, we use pho- 
tometry of the low-concentration, high-latitude (b=- 

6 The colour interval (mean colour and width) is compro- 
mise taken between incompleteness at faint magnitudes (and red 
colours) and the rising steepness of the MS in the CMD at blue 
colours. 
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47.7°) Galactic globular cluster Pal 12 from Martinez- 
Delgado et al. (2002; see Fig. E2 left). Pal 12 was ob- 
served through nominally the same filters and has a lim- 
iting magnitude similar to the CMa data. Taking E(B- 
V)sfd98 — 0.037 mag and a standard Galactic reddening 
law (see Sec. 12. 21 . we apply a single de- reddening value 
to the Pal 12 photometry. Note that we chose the Pal 
12 comparison for its matching filter and S/N; it does 
have a slightly different metallicity [[Fe/H](Pal 12)= - 
0.94dex, Harris 1996] than the old CMa MS ([Fe/H] ~ 
-0.7 dex, B04). As it is a very sparse cluster, we use the 
inner region (r< 1.5 arcmin), which reduces the field-star 
contamination. We then apply essentially the same pro- 
cedure to the Pal 12 photometry (Fig.^J right) as in the 
CMa analysis, except that we use a slightly wider colour 
range to have a well-defined MS. We obtain UMS.Pai 12 = 
0.16±0.02mag through data re-sampling. 

For the young over-density, we adopt two single- 
metallicity reference stellar populations spanning a wide 
age range as found for the young MS (Carraro et al. 
2006, B05, MD05). Specifically, we consider two metal- 
rich (Z=0.006) synthetic populations 7 , namely 1-13 Gyr 
and 0.1-13 Gyr, at the adopted D Q =7.5±lkpc, without 
binaries, with an adopted constant star formation rate 
and a zero distance range. 

6.2.1. Result - Depth of the old MS 

We compute the depth cr m in magnitudes as fol- 
lows: [0.49 2 - 0.16 2 ] 5 = 0.46 mag (± 0.05), where the 
(squared) values in parentheses are for observed MS 
width (see Sec. 16.1(1 and an estimate of the intrinsic width 
of a simple stellar population (see Sec. 16.21 above). The 
photometry error is negligible based on Fig. El As this 
is an upper limit, we have <7 m < 0.46 mag. Finally, us- 
ing equation we obtain cti os < 1.8 ± 0.2 ± 0.2 kpc, 
where the errors stem from the uncertainty estimate for 
cr m and the adopted distance uncertainty, respectively. 
Alternatively, FWHM los < 4 kpc. 

As a check, we use data given in MD05 (their Table 1). 
Subtracting their tabulated (random) errors, in quadra- 
ture, from the total (observed) MS width, we obtain tr m 
= 0.45 mag. Applying Eqn. El we have cti os = 1.8 kpc, 
which is in excellent agreement with the estimate in the 
present paper. 

6.2.2. Result - Depth of the young MS 

We now attempt a first estimate of the l.o.s. depth 
of the young MS, in a similar way to the old CMa MS 
analysis. As illustrated in Fig. ^3 we determine the B- 
band width of the young MS in the colour range 0.2 < 
(B-R)o < 0.5 mag. For the 0.1-13 Gyr and 1-13 Gyr syn- 
thetic populations, the estimated B-band width is 0.37± 
0.17mag and 0.51 ± 0.05 mag respectively. This suggests 
that there is a dependence on the age of the youngest 
synthetic stars, and a different metallicity may alter this 
further. The B-band width of the observed young MS is 
0.38±0.17mag. To estimate a preliminary upper limit on 
the l.o.s. depth, we consider the observed MS width plus 
its uncertaintly (i.e. 0.38 + 0.17mag), and subtract the 
width of the synthetic MS in quadrature. This gives cr m 
= 0.39 mag and 0.17 mag for the 0.1-13 Gyr and 1-13 Gyr 

7 Synthetic photometry was obtained using the IAC-STAR code 
(Aparicio & Gallart 2004). 



populations respectively. Using Eqn.[5]we obtain eri os < 
1.5 kpc and 0.6 kpc respectively. Considering that a bi- 
nary fraction and a metallicity spread have been omitted, 
which would decrease the estimated physical depth fur- 
ther, a strong implication is that the l.o.s. depth of the 
young over-density is eri os < 1.5 kpc. 

7. DISCUSSION 

7.1. Comparison with previous work 

We have obtained and analyzed a set of 59 CMDs 
towards the CMa over-density and overall we have di- 
rectly presented a representative sub-set of 18. Towards 
CMa, these CMDs reveal a prominent over-density of 
old and young MS stars at a defined distance range. 
In Fig. II If bottom-right), one can see that the number- 
density of old MS stars over the expected MW back- 
ground in the colour-magnitude plane is substantial, 
more than a factor of two for most (> 70%) of the old 
MS width. The density profile of old MS stars shown in 
Fig. El provides unambiguous evidence for a stellar over- 
density towards Canis Major. This validates the 2MASS 
red-giant-based detection of an overdensity in this direc- 
tion by Martin et al. (2004a). We find CMa to be sub- 
stantially more elongated in longitude than what Martin 
et al. found, in general agreement with B06 (see also 
Rocha-Pinto et al. 2006), whose work was also based on 
red giants from the 2MASS catalog. 

We also presented a first extended map that traces the 
surface density map of the young MS stars towards the 
CMa stellar over-density (Fig. El- While a kinematical 
link between the old and young MS populations is lack- 
ing, the map shows they are roughly co-spatial on the 
sky and that the younger of the two stellar populations 
is less spatially extended in both longitude and latitude. 

Fig. EI right) shows a comparison of the young and old 
MS profiles together with the red clump surface den- 
sity profile at b~-8° from B06 (estimated from their Fig. 
9). From Table 2 in B06 the density of red clump stars 
reaches a maximum near 1=244° which appears to be 
compatible with the possible peak in the young MS den- 
sity profile. No such peak is obvious in the old MS data, 
which are consistent with a near-flat profile in longitude 
across the entire survey width. Qualitatively, the old 
MS profile appears to be compatible with the result of 
Rocha-Pinto et al. (2006): a large, low-latitude stellar 
body whose surface density extends beyond 1> 270°. In 
this scheme, the 'original' CMa over-density (Martin et 
al. 2004a, B04, MD05) would be an outlying part of 
a larger over-density. Overall, it is noteworthy that we 
detect an over-density of both young and old MS stars 
in the same direction as the red-clump star-count pro- 
file given in B06 (their Fig. 9). Yet, it is unclear what 
actual relationship, if any, exists between them. Further 
study is required to place these issues on a firm statistical 
footing. 

We placed an upper limit of <7i os <1.8± 0.3 kpc 
(FWHMi os < 4 kpc) on the line-of-sight depth of the 
(old) CMa stellar over-density and found that it is highly 
elongated in projection (Al:Ab > 5:1). A comparison 
with the nearest known dwarf galaxy (Sgr), whose RR 
Lyrae star-based l.o.s. FWHM estimate is 5.3 kpc, cor- 
responding to eri os ~2.3kpc (deduced from (m-M)o data 
in Cserenjes et al. 2000) only suggests that the CMa 
depth of a few kilo-parcsecs estimated in this paper is 
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not necessarily incompatible with an origin as a dwarf 
galaxy. We also found that the depth of the young MS 
population, <7i os < 1.5 kpc, is apparently smaller than 
that of the older stellar population. The overall picture 
portrayed by the data in this paper is one of a young 
stellar population that is less extended, both in terms of 
its line-of-sight depth and angular size (both in longitude 
and latitude), than the older population. In Sec. 17. 21 bc- 
low, we scruntinize this briefly with regard to different 
interpretations of the CMa over-density. 

7.2. On the different interpretations of the CMa 
over-density 

As outlined in the introduction, there are three broad 
classes for the different interpretations of the CMa over- 
density, namely (a) it is a predominantly old (several 
Gigayear-old) dwarf galaxy that may be partially dis- 
rupted (Bellazzini et al. 2004, 2006; MD05); (b) it can 
be explained by a line-of-sight crossing a warped, but lo- 
cally axisymmetric (i.e. sub-structure-free) outer Galac- 
tic disk (Momany et al. 2004, 2006; Lopez-Corredoira 
2006); and (c) the over-densities in the young and old 
MS stars arise from out-of-plane MW spiral arms and 
the Local Arm respectively (Carraro et al. 2005, Moit- 
inho et al 2006). We discuss these in turn. 

(a) If the strongly elongated old over-density that we 
have mapped is the central part of a dwarf galaxy, could 
it still be partially bound? In a disrupting satellite, stars 
are lost from the gravitationally bound main body and 
carried in to tidal tails, one leading the satellite and one 
trailing it. If a satellite is on a near-circular, low-latitude 
co-planar orbit, as the CMa dwarf would have to be (e.g. 
Peharrubia et al. 2005), the predominant tidal forces are 
always in the direction towards the Galactic centre. This 
leads to perturbations over the entire orbital period, al- 
lowing stars to continually escape, unlike eccentric orbits, 
as evidenced by the modeled disruption of the globular 
cluster Pal 5 in Dehncn et al. (2004). But in such tidal 
disruption events (e.g. Piatek & Pryor 1995, Dehnen et 
al. 2004, Peharrubia et al. 2005), the still bound portion 
of the stars does not become highly flattened. The tidal 
debris, however, would be wrapped around the Galac- 
tic centre, as modelled by Peharrubia et al. (2005) and 
Martin et al. (2005). Based on these generic simula- 
tion results, we infer from the strong elongation of > 5:1 
in longitude that CMa has been undergoing recent tidal 
disruption and is still localized, but no longer gravita- 
tionally bound. 

Support for the dwarf galaxy hypothesis has come from 
the proper motion analysis in Dinescu et al. (2005), 
which is based on young MS stars in a small area 
(0.25 deg 2 ) near the suspected centre of the CMa over- 
density. It reveals a large (7a) deviation from the ex- 
pectation for stars belonging to the warped part of the 
Galactic disk, a motion that hts the pre-existing dynami- 
cal model (Peharrubia et al. 2005), that is itself tied only 
to known parts of the Monoceros stream. 8 This argu- 
ment of course hinges on the assumption that the young 
tracer stars are in fact associated with the predominantly 

8 We note that modeling of the progenitor of the Monoceros 
stream is discussed in Martinez-Delgado et al. (2005b), and shows 
a comparison of the N-body Monoceros model with the integrated 
orbit of CMa. 



older, larger and more complex CMa stellar over-density, 
an issue that is debated in Carraro et al. (2005). Simula- 
tions of dwarf galaxies in dynamical equilibrium suggest 
that young stars may well be kinematically colder than 
older populations (McConnachie et al. 2006). However, 
to firmly establish a link between the young and old CMa 
stellar populations or the absence of it, kinematical (and 
chemical) study of both would be essential. 

There is one aspect of the projected old and young 
MS over-density distributions that is not easily recon- 
ciled with the idea of both arising from a recently dis- 
rupted satellite: the large extent of the young popula- 
tion in longitude, given its more compact distribution in 
latitude. If the young stars were at the centre of the pu- 
tative precursor dwarf galaxy (as seen in e.g. Phoenix, 
Martinez-Delgado et al. 1999 and Fornax, Stetson et al. 
1998), then they should have been disrupted last, and 
hence the least 'stretched out' sub-population. Only yet 
more extensive imaging (in area coverage) can resolve 
this issue. 

(b) Could the old CMa stellar over-density simply be 
a consequence of viewing the warped outer stellar disk 
(Momany et al. 2004, 2006) nearly edge-on ? The 
2MASS red-giant star surface density map of Fig. 9 
(bottom-panel) in Momany et al. (2006) appears to 
be broadly compatible with the flat or possibly gradual 
density variation across our survey width from l=231°to 
258°. The angular density distribution of CMa that we 
have mapped appears, taken by itself, to be qualitatively 
consistent with a warped Galactic disk (Momany et al. 
2006). From their analysis (their Sec. 4.2), one would 
infer however that at D©=7.3kpc star-counts at b=+8° 
should match those at b^ -14°. This is not reflected in 
our data, as is evident from a comparison of the CMDs 
for b=+8° and -15° in our Fig. |SJ Rather, an additional 
stellar population is present at 1 = 230° to 260°, extend- 
ing up to 2 kpc below the Galactic mid-plane. 

More importantly, our line-of-sight depth estimate can 
shed new light on this matter. In Fig. [21 we compare 
the l.o. s. distribution predicted by the model in Lopez- 
Corredoira et al. (2002; LC02) and Yusifov (2004) to 
the one observed here. This figure shows (i) that the 
maximum reached by the old stellar over-density is sig- 
nificantly farther [A(m-M)~4.8 mag or ~ 10.5cr o id ms] 
from the maximum expected for a locally axisymmetric, 
warped and flared MW disk, based on LC02 (without 
spiral arms or other sub-structure). The distinctiveness 
of the old stellar over-density is still evident if we com- 
pare with the pulsar-based Yusifov (2004) model, which 
yields [A(m-M)~2.8 mag or ^6 0" o ldMs]- In the Yusi- 
fov (2004) model, the density profile has a more grad- 
ual drop-off with galacto-centric distance, and therefore 
reaches a maximum at larger helio-centric distances than 
the LC02 model. Fig.^^also shows (ii) that the observed 
CMa depth is much smaller than expected for a line-of- 
sight intercepting a warped disk. We conclude that the 
existing 'smooth' (i.e. locally axisymmetric) descriptions 
of the warped outer Galactic disk are unlikely explana- 
tions of the CMa over-density. As the ai os discrepancy 
(distance and width) seems generic, variations on the ex- 
isting warp models are also unlikely to match all data. 

Fig. 1151 reiterates this point by showing a colour- 
magnitude comparison of an observed CMa field, a con- 
trol field and a synthetic field for a warped and flared 
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MW stellar disk. The synthetic population is 4-10 Gyr- 
old with -0.4 > [Z/Zq] > -0.5 dex, and has been con- 
volved with the LC02 MW density profile plotted in 
Fig. 1141 It is in fair qualitative agreement with Fig. 7 in 
B06, who used a simple stellar population (47Tuc) as the 
template population. Our model CMD is not meant to 
be an exact simulation of the MW's contribution to CMa 
CMDs, but rather it serves to illustrate that 'smooth' 
(locally axi-symmetric) components of the warped outer 
Galactic disk alone do not account for the CMa over- 
density. 

In the model MW CMD, Fig. HU right), there is 
a relative increase in star-counts at [(B-R)~0.8 mag, 
19<B<17] that is not observed in the control field. It 
might be a signature of viewing the warped outer Galac- 
tic disk in projection, and differs substantially from the 
CMD feature known as the CMa stellar over-density, 
which is evident in Fig. ITBT left). 

Lastly, we note a possibly strong low-latitude (b=- 
3.4°) detection in Bragaglia et al. (2006) of the CMa 
over-density towards the Galactic anti-centre where the 
Galactic warp is largely absent. This could be related to 
a tidal tail of the (candidate) dwarf galaxy, and appears 
to disfavour a projection effect of the warped outer MW 
disk. 

(c) There are interpretations of the young MS pop- 
ulation that do not invoke a satellite origin, but in- 
stead infer that known types of MW sub-structure may 
have been detected (Carraro et al. 2005, Moitinho et al 
2006). In this scenario the young stellar over-density is 
a < 100 Myr-old spiral arm population and the old over- 
densities are merely a projection effect of looking along a 
nearby inter-spiral arm structure (the Local Arm). While 
we cannot provide firm evidence against the inter-spiral 
arm scenario, it is unclear whether it can be reconciled 
with the old stellar population, which is still detected 
at b^-15° in our survey, and is distributed over a larger 
area than the young stellar population. A comparison 
of the bluest stars in the empirical and synthetic CMDs 
in Fig. El suggests that the most recent star formation 
activity could indeed have occured less than ~ 100 Myr 
ago, but this depends on reddening and the choice of 
metallicity and distance. Young (<100 Myr-old) stars 
are not necessarily incompatible with the MD05 study 
in which the age used for the youngest MS stars is very 
uncertain, because of the same degeneracy between the 
distance, the stellar population (age, Z), and foreground 
reddening. The youngest stars could be 100 Myr-old if 
the over-density is farther way (e.g. ~ 10 kpc), as needed 
for the Peharrubia et al. (2005) model, which places 
the progenitor of the Monoceros stream at larger dis- 
tances too. A reliable distance estimate, based on RR 
Lyrae stars for example, as well as spectroscopic esti- 
mates and extinction information are necessary to solve 
this controversy unequivocably. We note that young (e.g. 
100 Myr) stars themselves are not strong evidence against 
a dwarf galaxy: For example, the last episode of star for- 
mation observed in the Phoenix dSph was ~ 100 Myr ago 
(Martinez-Delgado et al. 1999), and similarly in Sextans 
A (Dolphin et al. 2003), and in NGC 205 (Butler & 
Martinez-Delgado et al. 2005). 

There is also the issue of why the young stellar over- 
density can be spotted 10° (up to 1.3 kpc) below the 
Galactic mid-plane (b=0°). If we are observing 100 Myr- 



old members of a MW spiral arm or an inter-spiral arm 
structure, such stars at D Q =7.5kpc would have drifted 
away from their parent stellar associations by a negligible 
amount (few x lOOpc), leaving them absent at b~-10°. 
This suggests that the young stars may have unusual 
kinematics, as was found for previous detections of young 
(< 100 Myr) stars several kilo-parsecs from the mid-plane 
(e.g. Rodgers et al. 1981; Lance 1988); such young out- 
of-plane stars have been attributed by those authors to 
possible accretion events, based on stellar kinematics. 

8. SUMMARY AND CONCLUSIONS 

We have presented initial results from our imaging sur- 
vey of the Canis Major stellar over-density, providing 
a large, representative sub-set of colour-magnitude dia- 
grams from our survey area, a depth estimation for the 
young and old MS populations, and a first analysis of 
their surface density distribution. In particular, our key 
findings are these: 

• Using 'old' MS stars, we can delineate an over- 
density of MS stars elongated along galactic longi- 
tude at a distance of D Q ~7.5 kpc. It coincides with 
the over-density of 2MASS red-giants discovered in 
Martin et al. (2004a), but our surface density map- 
ping reveals it to be markedly more elongated than 
initially thought. Its projected aspect ratio is prob- 
ably > 5:1, which is consistent with the more recent 
2MASS analyses (B06; Rocha-Pinto et al. 2006). 
We also map the angular distribution of the much 
bluer MS stars (the 'young' MS). 

• The distributions of young and old stars are ap- 
proximately co-spatial in projection, but the young 
MS density profile is markedly more localized, with 
a shallow maximum near (1, b) ~ (240°, -7°). 

• We report the clear detection of young and old MS 
stars up to 1.3 kpc and 2 kpc respectively below the 
Galactic midplane at 1^240°. 

• We derive an upper limit on the line-of-sight depth 
of the (old) stellar over-density, by assuming that 
its stellar population is simple (as, e.g. globular 
cluster Pal 12), and that its MS width soley reflects 
a distance spread. We obtain ai os < 1.8 ± 0.3 kpc 
(or ~ 4 kpc, FWHMi os ) at the adopted D = 7.5 
± 1 kpc. The young MS stars are consistent with 
ctios < 1.5 kpc. 

We discussed these results in the context of three broad 
explanations put forth in the literature for the CMa over- 
density, (a) We infer from the strong elongation of the 
over-density in longitude and simulations in the litera- 
ture that, if it is a satellite galaxy on a near-circular, 
low-latitude orbit, it is unlikely to be still gravitation- 
ally bound at the present epoch; it would have to be a 
recently disrupted satellite, (b) The distance and line- 
of-sight depth of the over-density is in disagreement with 
all published 'smooth' (or locally axisymmetric) models 
of the Galactic warp. Those produce a density profile 
that is markedly more extended along the line-of-sight 
and reaches a maximum significantly closer. 

Lastly, of an out-of-plane spiral arm hypothesis for the 
young MS stars, we note that the presence of young out- 
of-plane stars is not uncommon in the MW, and such 
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stars in the literature have been attributed to possible 
accretion events, based on stellar kinematics. 

Without detailed modeling the data themselves are not 
yet sufficient to discriminate between an interpretation as 
sub-structure in the pre-existing warped outer Galactic 
disk or a disrupted satellite. 

It is a great pleasure to thank T. Erben for his very 
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This figure is given separately as a JPG file. 

Fig. 1. — Colour-magnitude diagram (CMD) for a field near the presumed centre of the CMa over-density, at (l,b) = (240.5°, -6.8°). It is 
labelled to show several different stellar components, namely the young and old CMa main sequences, and the red plume of late type thin 
and th ick d isk MS s tars . The approximate saturation level is marked (dashes). The CMD photometry has been de- reddened, as detailed 
in Sec. 12.21 See Sec. I4.ll for further explanations of the CMD content. 
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Fig. 2. — Sample completeness diagrams for the CMa photometry from artificial star tests as a function of the input artificial star 
magnitude. Coupled with Fig. [3] this is a diagnostic of the photometric quality. We have restricted the tests to six fields that provide 
a good coverage of the range of latitudes and longitudes of the CMa over-density. As these tests were performed using instrumental 
magnitudes, we applied an approximate zero-point correction to each bandpass to shift the magnitudes to the flux calibrated system. Also, 
as we use de-reddened CMDs in this paper, we applied an appropriate shift to the input magnitude scale of each completeness curve, 
denoted by Magnitude^. We refer to the magnitude shifted B- and R-band completeness curves by B' and R/ respectively in the panels. 
See Sec.|2|for further details. 
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These figures are given separately as a JPG file. 

Fig. 3. — Photometry diagnostic, based on artificial star tests. In each of the two groups of panels we provide magnitude scatter plots 
(upper) [recovered minus input magnitudes] as a function of the input magnitude, as well as the magnitude dependence of the error curve 
enclosing ±34% of the stars (lower panels). They indicate that 68% of injected stars at B'~22.6mag and R'~22mag in these fields have 
a photometry error typically below 0.05 mag, and smaller still at brighter magnitudes. For the sake of clarity, we restricted this diagnostic 
check to a sub-set of six fields spanning the survey area, with some at similar positions to check the local and global photometric quality. 
For the meaning of BJ N and RJ N , please see the meaning of MagnitudeJ N in Fig. |21 
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This figure is given separately as a JPG file. 

Fig. 4. — CMD of an example 'control field' at (l,b)=(239.9° ,-20°) from our survey (top) and its interpretation in terms of a synthetic 
MW model (bottom), which accounts only for the smooth and symmetric components (thin, thick disk, and halo). It shows that such fields 
comprise stars over a range of distances without prominent sub-structure along the linc-of-sight. There is an absence of stars at Bo < 14.5 
in the control field owing to detector saturation, and at Bo > 22.5 because of the magnitude limit. It can be seen in the Galactic model 
CMD that fainter MS stars at a given colour are farther away. The marked density of model stars at [Bo > 22, (B-R)o ~ 0.6-0.8 mag] in 
the synthetic CMD are ancient, metal-poor MS turn-off stars in the outer Galactic halo, but their number-density may be inaccurate. 
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Fig. 5. — CMD of example control fields (right) at b =(+f 5°, +8°, +5°) from our survey (top-to-bottom), together with the CMDs for b 
=(-15 ,-8°,-5°) (left). It illustrates that the structure of the CMDs for fields above and below the mid-plane differ over the full latitude 
range at 1=240° . 
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This figure is given separately as a JPG file. 

Fig. 6. — We plot a sequence of CMDs across the Galactic plane along a fixed longitude (1~24Q°), running from b = +15° to b = -20°. 
The text in each panel refe rs to the (l,b) position. The axes are labelled [B Q , (B-R)q] because the photometry was de-reddened in two 
different steps (see Sec 12.21 . Exceptions are the fields at b=+15° ,+8° ,+5° and -20°, whose correct axis labelling is [Bo, (B-R)o], but is 
excluded for the sake of clarity. The top-right panel shows photometry error bars (based on Fig. |3J at (B-R) =0, which are typical for 
the bulk of these CMD s. We also show the de-reddening vector arising from the additional (i.e. post-SFD*) de-reddening step, which 
is described in Sec. 12.21 The 80% completeness limit (dashes), is included in two panels, illustrating the faint limiting magnitude of the 
photometry. 
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This figure is given separately as a JPG file. 

Fig. 7. — As Fig. IH1 b ut for a sequence parallel to the Galactic plane at latitude (b~-8°), running from 1 ~ 234 to 1 ~ 258°. For other 
explanations, see Fig. 161 
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This figure is given separately as a JPG file. 

Fig. 8. — Illustrating the extraction boxes used for the young and old MS star count estimates. The box to the right of the old MS is a 
reference box of MW stars. The 80% completeness limit is marked. See Sec. l5.ll for more on the extraction boxes. 
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Fig. 9. — Mapping the CMa stellar overdensity of young and old CMa MS stars as a function of Galactic latitude (left) and longitude 
(right). Only fields with median E(B-V)sfd» below 0.30 mag are included. The open diamonds represent the density of old MS stars (as 
estimated in Sec. |5]and recorded in Table HI): the triangles represent the young MS stars. In addition, CMa density estimates from red 
clump data (Bellazzini et al. 2006) are shown in the right panel, based on their Fig. 9, together with an assigned 20% error. The maximum 
of the Bellazzini et al. values has been scaled to match our results. 
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This figure is given separately as a JPG file. 

Fig. 10. — Surface density map of the area of sky sampled, including high extinction fields (median E(B-V)sfd* > 0.30 mag; diamonds). 
For the sake of clarity, the WFI field-of-view is over-sized. We note that the pointings at 1>251° were originally designed to be control 
fields. The (relative) number-density of old CMa MS stars (bottom) and young MS stars (top) is intensity-coded. One sees an over-density 
of old MS stars that is most likely centred below the Galactic mid-plane (also see Fig.[SJ, and exhibits a pronounced elongation in Galactic 
longitude. The young MS population is probably also centred below the mid-plane, but is substantially more localized in the longitude 
direction. 
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Fig. 11. — Estimating the line-of-sight depth of CMa: We show the star-count distribution (left panels) and the associated CMD (right 
panels), for a reference star-field (top) and a field, (l,b) = (242.5°, -9.°) (bottom), near the presumed centre of the CMa over-density. 
The 'thickness' of the old MS, reflecting the l.o.s. extent is estimated in the colour interval bounded by vertical lines (dashes) in the 
bottom-right panel. As a guide to the photometric quality, the 80% completeness limit is also marked (horizontal and tilted dashes). In 
the bottom-left panel, the functional form from the top- left panel (short-dashes), together with a Gaussian function is fitted (long- dash es) 
to the histogram. The data in this panel are representative, and correspond to only one iteration of data re-sampling. See Sec. IQfor 
further explanations of the depth estimation. 
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Fig. 12. — Estimati ng t he MS width of a reference stellar population, the globular cluster Pal 12. We have applied a similar procedure 
to that shown in Fig. 1111 We extracted photometry from a narrow colour interval in the CMD (left-panel). We then fitted a Gaussian 
distribution (dashes, righ t-panel). The 1 a width data given in the right-panel is from one iteration of data re-sampling. For more 
information, see Sec. 16.21 
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Fig. 13. — Demonstrating the depth estimation metho d for the young MS population: We use same CMD (l,b = 242.5°, -9°) used to 
estimate the depth of the old CMa MS population (see Sec. l6.lt . We show the young empirical MS (top- left) and a corresponding histogram 
of star counts (top-right) from one iteration of data re-sampling, together with a fitted Gaussian distribution (dashes). The other pane ls 
are for the adopted synthetic stellar populations. Additional details on the depth estimation can be found in the panels and in Sec. 16.21 
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Fig. 14. — Illustrating that none of the 'smooth' (locally axisymmctric) MW models for (l,b) = (241.5° ,-7.5°) reach a maximum density 
near the distance modulus of the old CMa stars (solid). The area under each profile has been normalized to unity. The MW models 
are the warped, flared and inhomogeneity-free Galactic disk model from Lopez-Corredoira ct al. (2002) (short dashes) and the Yusifov 
(2004) model (long dashes). Model parameters are as given by those authors, except that we curtail the amplitude (z w ) of the warped 
Lopez-Corredoira et al. disk from diverging to infinity by setting z w (R>14kpc) = z w (R=14kpc), where R is galacto-centric distance. 
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This figure is given separately as a JPG file. 

Fig. 15. — Illustrating that the CMa stellar over-density at (l,b) = (240°,-8°) (left) is absent in the control field at b=+8° (middle), 
and is also absent in the (unreddened) Galaxy model CMD for b=-8° (right). The Galactic model has been constructed by convolvi ng a 
synthetic population (age =4-10Gyr, -0.4 >[Z/Zq]> -0.5dex), with the LC02 line-of-sight density profile of MW stars given in Fig. 1141 
Even though it has not been convolved with photometry errors, and contains no MW halo stars, a comparison of all three CMDs strongly 
disfavours the warped MW disk as a way to explain the old stellar over-density found in the CMa field. As we do not model the red plume 
of nearby (< 1.5 kpc) late-type MS stars, the red boundary is set to (B-R)=2mag. 
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a (l) Field ID; (2) and (3) Galactic longitude and latitude respectively; (4) Median differential extinction and standard 
deviation: a For this, data was taken from the Schlcgcl ct al. (1998) dust maps, modified using the correction from 
Bonifacio ct al. (2000; their Eqn. 1). (5) The additional colour dc-reddening, detailed in Sec. m\ (6) Num ber of stars 
in the CMD; (7) Number of CMa main sequence stars determined using the method described in Sec. nm (8) Number 
of young MS stars and the rms uncertainty, estimated inside in the young MS extraction box outlined in Fig. [8] There 
is no Noid MS or Nvoung MS estimate if E(B-V)sfd98* > 0.3 mag. We note that fields with negative Noid MS estimates 
are set to zero for the density profiles in Fig. [^1 and the surface density map in Fig. Ilul 



CMa stellar over-density 



TABLE 2 
Target information. Continued 
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a (l) Field ID; (2) and (3) Galactic longitude and latitude respectively; (4) Median differential extinction and 
standard deviation: a For this, data was taken from the Schlegel et al. (1998) dust maps, modified using the 
correction from Bonifacio et al. (2000; their Eqn. 1). (5) The additional colour de-reddening, detailed in Sec. 
(6) Number of stars in the CMD; (7) Number of CMa main sequence stars determined using the method described 
in Sec. I5.l1 (8) Number of young MS stars and the mis uncertainty, estimated inside in the young MS extraction 
box outlined in Fig. [8] There is no Noid MS or N Young MS estimate if E(B-V)sfd98* > 0.3mag. We note that 
fields with negative Noid MS estimates are set to zero for the density profiles in Fig. Inland the surface density map 
in Fig. EH 
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